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Graphene nanoribbons are the counterpart of carbon nanotubes in graphene-based nanoelec- 
tronics. We investigate the electronic properties of chemically modified ribbons by means of den- 
sity functional theory. We observe that chemical modifications of zigzag ribbons can break the 
spin degeneracy. This promotes the onset of a semiconducting-metal transition, or of an half- 
semiconducting state, with the two spin channels having a different bandgap, or of a spin-polarized 
half-semiconducting state — where the spins in the valence and conduction bands are oppositely po- 
larized. Edge functionalization of armchair ribbons gives electronic states a few eV away from the 
Fermi level, and does not significantly affect their bandgap. N and B produce different effects, de- 
pending on the position of the substitutional site. In particular, edge substitutions at low density do 
not significantly alter the bandgap, while bulk substitution promotes the onset of semiconducting- 
metal transitions. Pyridine-like defects induce a semiconducting-metal transition. 



I. INTRODUCTION 

Graphene is the latest carbon allotrope to be discov- 
ered [l|, |2|, |j, |5j , and it is now at the center of a sig- 
nificant experimental and theoretical research effort. In 
particular, near-ballistic transport at room temperature 
and high carrier^mobilities (between 3000 and 100000 
make it a potential material for na- 
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applications. 

It is now possible to produce graphene samples with ar- 
eas exceeding thousands of square microns by means of 
micro-mechanical cleavage of graphite, and much larger 
by " expitaxial " growth on SiC , @, E, [EE H, [H, [H, E3] ■ 
An ongoing effort is being devoted to large scale produc- 
tion a nd g rowth on different substrates of choice [IsL EE 

Emm. 

Graphene can be readily identified in terms of num- 
ber and orientation of the layers by means of elastic 
and inelastic light scattering, such as Raman |23l. HE 
HE HE HI and Rayleigh spectroscopies [H, H3T Ra- 
man spectroscopy also allows monitoring of doping and 
defects (23l . I25I l30l. [3l| . Once identified, graphene lay- 
ers can be processed into nanoribbons by lithography 
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As for carbon nanotubes (CNTs) [ISS], electron con- 
finement modifies the electronic structure of graphene, 
when this is cut into nanoribbons (GNRs) [3HL l3a . I37L [H, 
I39I HE EH . When GNRs are cut from a single graphene 
layer, their edges could in general consist of a combina- 
tion of regions having an armchair or a zigzag geome- 
try [IE HE H, 53 • If a ribbon is uniquely limited by 
one of these edges, it is defined either as an armchair 
GNR (AGNR) or as a zigzag GNR (ZGNR) (See, e.g., 

Fig.rr]) [IE ESS!. 

GNRs are the counterpart of nantubes in graphene 
nanoelectronics. Indeed, the confinement of the elec- 
tronic wavefunctions and the presence of the edges open 
a bandgap, making them suitable for the realization of 
devices. Due to their potential technological applica- 
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ticular attention to the factors determining the presence 
and the size of the gap. 

CNTs can be metallic or semiconducting, depending on 
their chirality. This could lead to a fully carbon based 
electronics, where semiconducting tubes are used as 
channels and metallic ones as interconnects [HI. IriE l66j|. 
However, the present lack of control on the chirality pre- 
vents to engineer their electronic properties on-demand, 
and is a major barrier to industrial implementation. 

GNRs could combine the exceptional properties of 
graphene, and the possibility of being manufactured by 
means of industry-amenable large scale top down planar 
technologies, which would uniquely define their chirality. 

Introducing impurities and functional groups can be an 
effective way to control electronic properties of GNRs. 
On the other hand, covalently bonded impurities are 
likely to be a by-product of the production processes. 

Here, we investigate the effect of chemical disorder on 
the band structure of ZGNRs and AGNRs. We first study 
edge functionalization of ZGNRs with a set of different 
radicals, elucidating the role played (i) by the ribbons 
width, (ii) by the concentration of the functional groups 
along the ribbon's edges, and (iii) by the presence of 
one- and two-sided edge functionalization. We find that 
edge functionalization breaks the spin degeneracy, lead- 
ing to four possible outcomes. In the first case, the rib- 
bon maintains its semiconducting nature, and the top of 
the valence band and the bottom of the conduction band 
belong to the same spin channel. We refer to this as 
"spin-selective half-semiconductivity" In the second 
case, the ribbon is still semiconducting, but the top of 
the valence band and the bottom of the conduction band 
belong to opposite spin channels. We refer to this as as 
"spin-polarized half-semiconductivity" p71 |. In the third 
and the fourth cases, the band gap of one or both spin 
channels closes. In the latter case, the ribbon undergoes 
a semiconductor-metal transition, while in the former, it 
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behaves as a half-metal [50, 68]. Edge functionalization 
on AGNRs at the moderately low densities we studied 
has negligible effect on the gap. 

We then consider chemical doping on H-terminated 
AGNR and ZGNR. AGNRs are always semiconducting, 
so, in principle, B, N, or atomic substitutions can re- 
sult in some form of electrical doping. However, given the 
particular nature of the ribbons, it is not clear a-priori 
whether chemical doping of AGNRs would result in the 
formation of defect states within the gap, like in carbon 
nanotubes (69l . [toI . l7ll . E2] and bulk semiconductors, or if 
it would rather cause a shift of the Fermi level similarly 
to what observed in graphene [H, [zE Iz3|- Thus, we in- 
vestigate in detail the effects of atomic substitutions on 
the electronic structure of AGNRs. We find that upon B 
and N edge-substitution, impurity levels appear far away 
from the band gap, slightly renormalizing its width. On 
the contrary, N and B substitution in the bulk of the rib- 
bons result in impurity levels near the gap, and in a shift 
of the Fermi energy within one of the continuum levels 
of the pristine ribbon. Atomic substitution in edge func- 
tionalized ZGNR and pyridine-like impurities in AGNR 
induce semiconductor-metal transitions for high impurity 
densities, but do not give impurity levels close to the top 
or bottom of the gap, unlike the case of nanotube doping. 



II. BACKGROUND 
A. Band structure of GNRs 

GNRs can be thought as single wall CNTs cut along 
a line parallel to their axis and then unfolded into a 
planar geometry. Since on a graphene plane the zigzag 
and armchair directions are orthogonal, this procedure 
transforms an armchair CNT into a ZGNR and a zigzag 
CNT into an AGNR. However, this correspondence be- 
tween the geometry of CNTs and GNRs is not entirely 
reflected in their electronic properties. Indeed, in con- 
trast to CNTs, where armchair tubes are always metallic 
and zigzag could be metallic or semiconducting depend- 
ing on the chiral angle ||75j, the earliest theoretical stud- 
ies of the electronic structure predicted all ZGNR to be 
metallic, while AGNR were expected to be semiconduct- 
ing [H, ES E3, El, ES E3|- The metallic character of 
ZGNR was attributed to the presence of a high density 
of edge states at the Fermi energy [3f| E3, El, 13] . 

More recently, spin-polarized DFT calculations have 
found that the ground state of ZGNR has a anti- 
ferromagnetic configuration ( AF) , where electronic states 
with opposite spin are highly localized at the two GNR 
edges [36|,EilE3,El,E§] and are responsible for the open- 
ing of a gap. It was shown that the magnetic instability 
energy of ZGNRs, i. e. the energy difference between 
the ferromagnetic and the non magnetic ground state, in- 
creases from ~ 0.25 eV/unit cell up to ~ 0.37 eV/unit cell 
when N increases from 5 to 30, and then stabilizes [4o| . 
The energy difference between the ferromagnetic and the 



AF ground state, however, are much smaller, with the 
AF state being more stable than the ferromagnetic by 
4.0, 1.8, 0.4meV/unit cell when N is 8, 16, 32 PS]. Since 
UbT at room temperature corresponds to ~25 meV, this 
indicates that ZGNRs at room temperature are stable 
in a magnetic state. As temperature decreases the anti- 
ferromagnetic state becomes favored over the ferromag- 
netic one. Thus DFT predicts both AGNRs and ZGNRs 
to always have a direct nonzero band gap at least at low 
temperatures [H, E3, IHoll . In CNTs the gap scales in- 
versely with diameter [7(j. Similarly in GNRs the gap 
scales inversely with the width [H, Eo, EH, EH, E3|, clos- 
ing, as expected, for infinite graphene. 

The electronic structure of GNRs was also studied us- 
ing GW calculations [HSUl- As for CNTs Jz3|, the better 
description of the electron-electron interaction provided 
by GW increases the computed GNRs electronic gap. 
Furthermore, calculations based on the Bethe-Salpeter 
equation proved that, due to the presence of excitons, in 
GNRs the optical g ap is significantly smaller than the 
single particle one [H, EH- However, apart from the 
gap renormalization, all the main results obtained within 
DFT, including the AF ground-state of ZGNRs and the 
width dependence of the bandgap, are confirmed by the 
computationally more expensive GW [H, E3| ■ 

The presence of a gap inversely proportional to the rib- 
bon's width was experimentally observed, for ribbons in 
the 10 to 500nm range, by means of temperature depen- 
dent conductance (G) measurements (1,[1, [ll]. The gaps 
(Eg) measured in these experiments are larger than those 
predicted by theory, and their width scales as E g =A(W- 
W*)^ 1 , with W being the ribbon's width and A and 
W* fit constants [1]. In the same samples, the conduc- 
tivity G was found to scale as G=B(W-Wo), with B 
and Wo fit constants. For T=1.6K, it was found that 
W*^Wo=16nm @. This suggests that structural dis- 
order at the edges caused by etching, or inaccuracies in 
the width determination due to over-etching underneath 
the mask, can result in an "effective transport width" 
smaller than the ribbon's nominal width. Within this 
picture, Wo and W* can be considered as the reduc- 
tion of the effective width with respect to the nominal 
width Q. Based on these results, some authors have also 
proposed that the size of the GNRs bandgap can be ex- 
plained not only by confinement, but also by Coulomb 
blockade originated by the roughness at the edges of the 
ribbons 14411 . This would also determine the value of 
W^*~Wo [4J|. Furthermore, it was also found that the 
minimum conductivity as a function of the gate changes 
with the environment or the GNR thermal history [§]. 

Finally, the electronic structure of mono-atomic step 
edges of graphite was also probed by scanning tunnelling 
microscopy (STM) and spectroscopy (STS). This showed 
a large electronic density of states (DOS) at the Fermi 
energy (Ep) of zigzag steps, while a bangap was ob- 
served on armchair edges (42, EE EE, These re- 
sults seem to corroborate the simple tight-binding models 
predicting an accumulation of electronic states at Ep in 



ZGNRs [36j, [42, |43j, [4fl|, and to contradict the presence 
of the AF bandgap predicted by DFT and GW. How- 
ever, it is important to recall that such experiments were 
performed on monoatomic steps of bulk graphite, which 
correspond to a semi-infinite ribbon. Since the onset of 
an antiferromagnetic state requires the presence of wave- 
function with opposite spins localized on the two oppo- 
site edges of a ribbon, such state cannot be present in a 
semi-infinite system with just a single edge. 



B. Spin properties 

Ref. [50] showed, using DFT, that applying an electric 
field parallel to the ribbon's plane results in an opposite 
local gating of the spin states on the two edges of the 
ZGNR. The required critical field was estimated to be 
3.0V/W, WinA [53]. Indeed, the in-plane field lifts 
the spin degeneracy by reducing the band gap for one 
spin channel, while the other experiences a gap widen- 
ing. This promotes the onset of half-metallicity [68, 80] or 
induces a spin-selective semiconducting behavior, where 
the band gap of one spin is bigger than the other. Which 
of these states is established depends on the intensity of 
the applied field [50l . IBR IBtl and on the ribbon's length 
[56l. I57I. I5S| . These phenomena may have important ap- 
plications for the fabrication of spin filtering devices|8ll|. 

Half-metals are materials in which the current can be 
completely spin polarized as a result of the coexistence 
of a metallic nature for electrons with one spin orien- 
tation and an insulating nature for electrons with the 
other spin [H3, [HS] - Some ferromagnetic metals were pre- 
dicted to have this behavior [68| and were first observed 
experimentally in a manganese perovskite (sjj. Simi- 
larly, there can be half semiconducting states in which 
the electron spins in the valence and conduction bands 
are oppositely polarized [67]. Since the development of 
efficient devices for spin injection, transport, manipula- 
tion and detection is a fundamental requirement for the 
incorporation of spin into existing semiconductor tech- 
nology and to achieve spin based quantum computing 
and spintronics [HI], substantial efforts were made in to 
find half- metallic materials [fill, It is thus of great 
interest to identify such states in GNRs. 



C. Chemical modifications 

Chemical functionalization and substitutional doping 
have been investigated for many years in nanotubes [69 . 
[lE ED, E2, H3, H3, H3l, with the aim of tailoring their 
properties for sensing, transport, chemical and opti- 
cal applications or to incorporate them into polymers 
[H, [83, [sl, [H It is thus natural to do a similar 

investigation for GNRs. 

Various theoretical studies of the influence of edge 
geometry and chemical modifications on the transport 
properties of ribbon junctions have been reported so 
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Figure 1: Scheme of H terminated (a) zigzag, and (b) arm- 
chair nanoribbons. Periodic boundary conditions are assumed 
in the x direction. N and M are the numbers of "colums" and 
"rows" of atoms used to label the MxN ribbon. 



far [59j, [6CJ, |9l|, |92, |93j, |9J] . Furthermore, several works 
considered the reactivity of the ribbon edges and the ef- 
fects of geometrical and chemical modifications on their 
electronic properties [H, HI El, ED, M, El El El EE M, 
[99l . llOOL llOll . Il02l | . In particular, edge functionalization 
with oxygen-contai ning groups [55j and edge doping or 
edge imperfections [l03l | were shown to significantly lower 
the electric field required to induce half-metallicity in 
ZGNRs [55j], to break the spin symme try a nd to promote 
the possibility of net spin injection [l03j . B substitu- 
tional doping (621 . [63| , either on the edge or in the bulk of 
GNRs, was reported to cause metal-semiconductor tran- 
sitions in the ferromagnetic and non-magnetic ZGNRs 
states [filEo!]. The effect of B on the AF ground state 
has not yet been studied. 

Also of in teres t is the effect of edge radical functional- 
ization [99l . Il0lj |. Carbon atoms at the edges of ZGNRs 
offer superior chemical reactivity for the attachment of 
chemical groups than those on A GNR s edges, in the rib- 
bon's bulk, or on CNTs edges fool . 1 1 5| | . Indeed , as i n th e 
case of herringbone and bamboo-shaped CNTs |l06l . ll07| . 
the reactivity of GNRs is due to the presence of open- 
ended graphene sheets, which, in the case of ZGNRs, are 
particularly active thanks to the presence of electronic 
states highly localized at the edge carbon atoms [4^.[99|. 

We note that transport and Raman measurements 
show that as-prepared graphene samples and devices 
reach the charge neutrality point for gate voltages dif- 
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ferent from zero (H, ED, III, III, S Efll. Most times 
this corresponds to a shift of the Fermi level towards the 
valence band (indicative of p-doping) [13, 0, Il3|, but 
sometimes also towards the conduction band (equivalent 
to n-doping) (7£|. In particular, refs. (73l. l74| seem to in- 
dicate a bigger intrinsic p-doping for smaller ribbon size. 
Whether this behavior is a consequence of chemical dop- 
ing due to edge functionalisation during lithography, or 
of a rearrangement in the electronic structure due to the 
presence of defect states away from the Fermi energy, 
or, as for large grap hene samples, due to the presence of 
adsorbates [23, 108| needs to be further investigated. 



III. METHODOLOGY 



We perform spin polarized ab initio calc ulations with 
the CASTEP plane wave DFT code [l09| on hydrogen 
terminated ribbons. We use the Per dew-Burke-Ernzerhof 
gradient corr ected functional |llOl | and ultrasoft pseu- 
dopotentials [lllj with cut-off energies of 400eV. Rib- 
bons are built from the perfect graphite geometry with 
an initial C-C distance of 1.42 A. Geometry optimiza- 
tion is performed for all the structures with th e Bro yden- 
Fletcher-Goldfarb-Shanno (BFGS) algorithm [HI] until 
all the forces are smaller than 0.04 eV/A and the stress 
in the periodic direction is less than 0.05 GPa. 

A 4-atom unit cell is used to build our ZGNRs as shown 
in Fig.QJa) and the same cell rotated 90°is used to build 
the AGNRs of Fig. Q](b). The brillouin zones (BZ) of 
the GNR unit cells are sampled by Monkhorst-Pack [ll3| 
grids of the form Pxlxl with P such that the maxi- 
mum spacing between fc-points in the periodic direction 
x (Fig. [J) is < 0.1 A -1 . We find that, to simulate iso- 
lated ribbons, the in plane and perpendicular distances 
between ribbons in adjacent super-cells have to be larger 
than 5.5 A and 6.5 A, respectively. The band structure 
is calculated with an eigenvalue tolerance of 0.001 meV. 

We classify the ribbons with a MxN convention, 
where M and N are the number of rows and columns 
across the GNR width and length, respectively, as shown 
in Fig. [TJ Throughout this work, we refer to edge func- 
tionalized ribbons with the MxN index of the clean 
GNR preceded by the chemical symbol of the functional- 
izing radical and followed by the character Z or A denot- 
ing the zigzag or armchair chirality. Thus, e.g., NH 2 on 
the edge of a 2x5 zigzag ribbon is denoted as NH 2 -2x5Z. 

Atom substituted ribbons are labeled by the symbol 
of the substitutional atom followed by the word bulk or 
edge, to denote bulk or edge substitution, and the MxN 
index of the ribbon. For example N bulk substitution in 
the 6x5 AGNR is denoted as Nbulk-6x5A. 



IV. RESULTS 
A. Edge functionalisation: single edge 

Due to the presence of edge states in ZGNRs and 
their absence in AGNRs, the effects of single-edge radi- 
cal functionalization strongly differ. Indeed, as shown in 
the following and in Figs. EJa-f), [UJa), functionalization 
of ZGNRs lifts the spin degeneracy and in some cases 
also promotes a semiconductor-metal transition, while in 
AGNR the largest effect is the presence of impurity levels 
deep in the valence band that hardly modify the gap. 

We start by investigating five radical groups on the 
6x2 ZGNR. The top insets in Fig. Efa-f) schematize 
functionalisation with (a) NH 2 , (b) OH, (c) COOH, (d,e) 
N0 2 and (f) O. In general, geometry optimization shows 
that edge functionalization does not alter significantly 
the GNRs structure. For these radicals, we find the fol- 
lowing optimized geometries: the NH2 sits 1.386 A away 
from the edge carbon, the H-N distances are 1.021 A 
and the H-N-H plane is tilted ~20° with respect to the 
plane of the ribbon. The OH remains on the graphene 
plane with C-0 and O-H distances of 1.367 A and 0.977 
A respectively and a C-O-H angle of ^109°. The plane 
defined by the O-C-0 of the COOH has an angle of -53° 
with respect to the plane of the ribbon and the C-C, 
C-Ol, C-02 and 02-H distances are 1.487 A, 1.229 A, 
1.365 A and 0.985 A respectively. For NO2 we find two 
equilibrium geometries: one energetically stable (S) and 
one metastable (MS), separated by 0.407 eV. The S-N0 2 
has C-N and N-0 distances of 1.466 A and 1.257 A, and 
the O-N-0 plane is tilted ^60° with respect to the cross 
section of the ribbon. The MS-NO2 sits in a configura- 
tion similar to the NH2, with C-N and N-0 distances of 
1.466 A and 1.262 A. Finally, for the functionalization 
with a single O, this sits in the ribbon's plane in front of 
the edge C with C-0 distance 1.245 A. 

Functionalization strongly affects the electronic struc- 
ture of ZGNRs. This is shown in Fig. [2](a-f ) , where the 
density of states (DOS) of a set of functionalized rib- 
bons is compared with the corresponding DOS of the 
clean ribbon. Single-side radical edge functionalization 
of a ZGNR (Figs. I2I4H has an effect on the ribbon's spin 
density and band structure similar to that of an electric 
field parallel to the ribbon's plane[H3, HH, [H, Thus, 
the presence of functional groups, while preserving an 
AF ground state [39L [5l| , lifts the spin degeneracy of the 
clean ribbons [10l|. For the radical density we use, this 
results in the onset of spin-selective semiconductivity in 
NH 2 , OH, COOH and S NO2 functionalizations, and of 
a semiconductor-metal transition for O and MS NO2. 

The lifted spin degeneracy implies that in all cases the 
band gaps of the two spin channels, which we call a-spin 
and /3-spin, are different (Fig. |4] a- red/gray solid line 
and /3-blue /black solid line) . The change of gap with re- 
spect to the original H terminated GNR fingerprints the 
radical used to functionalize the ribbon. This is shown in 
Fig.[2ja-f), which plots the DOS for each of the five rad- 
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Figure 2: (Color online) Spin density of states of a (red/gray solid line) and /3 (blue/black solid line) for functionalized ZGNRs 
in comparison with clean ribbon (black dotted line), (a-f) 6x2Z (i. e. M = 6, N = 2 ZGNR) functionalized with (a) NH2, (b) 
OH, (c) COOH, (d) Stable-N02, (e) Meta Stable-NCh and (f) O radicals in one of the ribbon edge producing breakage of the 
spin degeneracy (a-f) spin-selective band gap(a-d) or semiconductor- metal transition (e,f). (g) Oxygen edge substitutional atom 
favoring a semiconductor-metal transition, (h) 6x5 zigzag ribbon with a nitrogen bulk substitution in the center of the ribbon, 
(i) Nitrogen bulk substituted ribbon shown in (h) with an extra single NH2 at the edge of the opposite carbon sublattice, with 
the maximum distance between N and NH2 in the periodic direction presenting semiconductor-metal transition 



icals. A pristine 6x2Z ribbon has a 0.58 eV direct band 
gap at k = 0.57r/c, c being the size of the ribbon's unit 
cell along its periodic direction. The different function- 
alizations change the band gap as follows. In presence of 
NH2 the gap remains direct and at k = 0.5n/c for both 
spins, but shrinks to 0.42 eV and 0.39 eV for the a and 
(5 spin channels, respectively. A similar effect is obtained 
by means of OH functionalisation, with final gaps of 0.57 
eV and 0.49 eV. In the case of COOH, the direct gap at k 
= 0.5tt/c increases up to 0.62 eV, but an indirect one of 
0.45 eV, with extremes at k = in the conduction band 
and k = 0.5n/c in the valence band, appears. Similarly, 
S-NO2 gives a direct gap of 0.66 eV and a smaller indirect 
gap of 0.28 eV with the extremes located at the same po- 
sition as the indirect gap due to COOH. MS-NO2 and O 
have several extra levels in the valence band and a zero 
gap for both spin channels. This makes them metallic. 

The semiconductor-metal transition we predict for O 
functionalized ZGNRs is an intriguing. Ribbons are of- 
ten produced by oxygen plasma etching (78|, suggesting 
that the same technique used to cut ribbons, could also 
influence their electronic properties. We thus perform 



a calculation for a lower edge functionalization density. 
We find that for a 6x6Z ribbon the semiconductor- metal 
transition does not happen, but we have a gap of 0.82 eV 
and 0.41 eV for the a and j3 spins, respectively (+0.40% 
and -30% with respect to pristine ribbon). This implies 
that the semiconductor/metallic behavior of the ribbon 
could be controlled by the degree of oxidation of a single 
zigzag edge, provided this does not revert to OH, which, 
as discussed before, would not significantly affect the gap. 

We note that, while NO2 and NH2 on graphene act as 
strong acceptor and donors [l08l . lll4| . the same molecules 
at the edges of the ZGNRs modify the electronic struc- 
ture in the gap vicinity, but do not act as dopants. 

Fig. [3] shows the effect of functionalization on the spin 
density and band structure of both single- and double- 
side functionalized ZGNRs. Black dotted lines represent 
the bands of the 2x1(4x1) clean ZGNRs folded into the 
2x2(4x8) reciprocal unit cell. The band structure of the 
clean ribbons is spin degenerated [H, HH, and the spin 
density is distributed in an AF arrangement from the 
edges to the center, as in Refs.[39l. IHlj. Fig[3|a,d) indi- 
cates that single-edge NH2 functionalization introduces 



6 




1.5 

0.75 
| 

Eu 

LU 

-0.75 



-1.5 



if 


SD 

».•••.•.< 


V p — 

clean ■■■■ 


> • 


p • I 




LUMO 




»•••••••' 

••••••••< 

. m 



k(2n/c) 0-5 



b) o . 




1.5 

0.75 
I 

5i o 

LU 





s 1 . 




i, • .».•.■ 


if 


>•••••••< 




> • i 


'-J 


H<jMO 


/\ 


» • 4 


n 

u 


LUfrlO . 
'••.••••« 







k(2itfc) 0.3 



c) 



1 5 



0.75 



-075 



-1 5 




LUMO 

••• 
••• ••• 



0.5 



0.75 



tij 
-0.75 



-1.5 





SD 
•,*,•, 

■»•••••»*••••••••• 

••••••••••*•*•*•* 

• '•••*•'•'•■•*•"• 


m 

■ • • • • *.'.< 






,•.•.•,« • 

• ••••• • 


LUMO 



k(27r/c) 



0.5 




0.75 
? 

LU 

-0.75 





SD • 
> . ■ i • > t B * a * 

.... 




HOMO • 
■ • m • 

« • • ■ • 






\m i 

»•••••••< 



k(2l!/C) 5 



f) 5 




k(2n/c) °' 5 



Figure 3: (Color online) Band structure of the NH2 2x2Z (a- 
c) and 4x8Z (d-f) single-side (a,d) and double-side (b,c,e,f) 
functionalization. Single side functionalization (a,d) breaks 
spin degeneracy, whereas this is nearly recovered for double- 
sided antisymmetric (b,e) and symmetric (c,f) functionaliza- 
tion. O-ifi) spins are represented by red/gray(blue/black) solid 
lines. The folded band structures of the pristine ZGNRs are 
plotted with black dotted lines. On the right of each BS are 
the corresponding spin density maps (/3-blue and a-red) of 
the total spin density (SD), spin density of the HOMO and 
LUMO bands. The image shows the relaxed structure in xy 
projection and a close up of the edges in the yz projection. 



extra charges that alter the spin density. Due to the 
zigzag geometry, C atoms sitting on the outmost position 
of opposite edges always belong to different sublattices. 
As a consequence, single edge functionalization affects 
prevalently the spin channel localized on the same edge 
where the functional group is attached. This is shown in 
Fig. [3ja,d). Indeed, after functionalization, the two sub- 
lattices are visible in the total spin density (SD) maps 
with the characteristic AF spin arrangement. However, 
NH 2 alters the SD, attracting a spin regions in its vicin- 
ity, while the (3 spin density is substantially unchanged. 

The perturbation of the SD is reflected in the lifting of 
the spin degeneracy illustrated in the DOS plot of Fig. [21 
and in the band structure of Fig. [3]Ja,d). As the highest 
occupied molecular orbital (HOMO) and lowest unoccu- 
pied molecular orbital (LUMP) are highly localized on 
the edge of the ZGNRs [H,[3§|, the influence of edge func- 
tionalization is predominantly on these orbitals. This is 
shown by the maps of the total spin density and of the 
HOMO and LUMO spin denisity in Fig. EKd). For the 
valence band orbitals of the 2x2 and 4x8 ribbons, the 
a spin (red/gray) hosting the NH2 presents a stronger 
modification than the (3 spin (blue/black) orbital in the 
opposite edge. Looking to the corresponding band, the a 
HOMO is more modified than the (3 HOMO with respect 
to the original HOMO (black dotted line). The breakage 
of the spin degeneracy at the top of the HOMO band is 
0.21 eV and 0.13 eV for the NH 2 -2x2Z and NH 2 -4x8Z 
ribbons respectively (note the color matching between 
spin density maps and spin bands in all Fig. [3|) . We ex- 
pect similar effects for any N, M combination and for any 
radical attached to only one ribbon edge. 

We finally consider edge radical functionalization on 
AGNRs. AGNRs do not have edge states near the 
band gap, so the effect is expected to be weaker than 
in ZGNRs. Fig. [5f a ) shows NH2 functionalization on the 
edge of the 12x3 AGNR. Here an impurity level is found 
~1.5 eV below the edge of the valence band, the gap in- 
creases by +0.02 eV (~4%), with no lifting of the spin 
degeneracy. Thus, for AGNR, NH2 edge functionaliza- 
tion does not significantly affect the electronic proper- 
ties, since the impurity level is located deep in the va- 
lence band and no other major change happens to the 
DOS. We expect that functionalization with other chem- 
ical groups will give similar results. Nevertheless, finding 
a chemical group able to introduce impurity levels in the 
gap will be of great interest for electronic applications, 
so further investigations are desirable. 



B. Edge functionalisation: double edge 

The effects of double edge functionalization are shown 
in Fig. [3](b,c,e,f). Interestingly, when another radical of 
the same species is located in the opposite edge of the rib- 
bon, the spin degeneracy is nearly restored . Fig. [3](b,e) 
show the 4x8Z and 2x2Z ribbons with two NH2 impu- 
rities on opposites edges in an asymmetric arrangement. 
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N M 

Figure 4: (Color online) Change of spin band gap (/3 channel-blue dotted line, a channel-red solid line) as a function of (a) edge 
density -N- and (b) width of ribbon -M- for different set of edge functionalized ZGNRs. 



The breakage of the degeneracy between a and (3 HOMO 
is only 1 meV and 0.4 meV for the double side asymmet- 
ric functionalized 2x2Z and 4x8Z ribbons, respectively. 
Fig. [3tc,f) show the same ribbons with two NH 2 impu- 
rities sitting on opposite edges in a symmetric arrange- 
ment. In this case the breakage of the spin degeneracy 
between a and (3 HOMO bands is 11 meV and 3.6 meV 
for the 2x2Z and 4x8Z ribbons respectively. 

The differences between the symmetric and antisym- 
metric geometries indicate the different stability of these 
configurations. Indeed, the total energy of the antisym- 
metric ribbon is lower than that of the symmetric ribbon 
by 42 meV and 34 meV for the 2x2Z and 4x8Z ribbons. 

To restore spin degeneracy in 2-sided functionalized 
ribbons, it is necessary to have the same impurity den- 
sity on both edges, and the maximum separation between 
consecutive radicals of opposite edges has to be achieved. 
Finally, it is important to stress that any geometry distor- 
tion or irregularity will result in a stronger magnetization 
of one of the edges, thus lifting the spin degeneracy. 



C. Edge functionalisation: density of radicals 
dependence 

We now consider the effects of the spacing between 
radicals, which is proportional to the index N, on the 
electronic structure of one-sided functionalized ZGNRs. 
Fig. |UJa) plots, for a set of NH 2 and N0 2 functionalized 
ribbons, the change in the band gap for each spin as a 
function of TV. A high level of NH2 functionalization of 
the M = 2 ribbon can shrink the band gap of the a spin 
up to ~60%, while leaving the (3 gap almost unchanged. 

A similar behavior is observed for the MS-N02-2xiVZ 
ribbon. In this case, increasing the functionalization den- 
sities narrows the a gap up to almost its complete clo- 
sure, while the {3 gap has a maximum change of <~-20%, 



almost independent of doping. These results, obtained 
for very thin ribbons (M = 2), seem to suggest a corre- 
lation between the breakage of the spin degeneracy and 
the edge-functionalization density, the former slowly de- 
creasing for increasing N. However, in wider ribbons 
(e.g. M = 4, 6), the band gap shows a much weaker 
dependence on N. In particular, for the NH2-4xiV rib- 
bon the changes in the a and (3 gaps oscillate between 
20-35% and 3-18% respectively when N spans between 2 
and 8. Almost identical gap variations are observed for 
NH 2 -6xiVZ for N between 2 and 5. 



D. Edge functionalisation: width dependence 

Fig. Hfb) plots the variation of the a and f3 bandgap, 
relative to the size of the bangap in the clean ribbon, 
in single-side functionalized ribbons as a function of M . 
Fig.[4]suggest that the modulation of the gap due to func- 
tionalization depends very on M, although some strong 
fluctuations may still occur for the thinnest ribbons (i.e. 
M = 2). The points for M = 15 in Fig.gfb) confirm this 
behavior, showing that the variation of the a-spin gap 
due to functionalization in the S-N02-Afx2Z ribbons is 
-32% to -55% for M changing from 2 to 6. A further in- 
crease of M from 6 to 15 (ribbon width increasing from 
1.3 nm to 3.8 nm) leaves the functionalization contribu- 
tion almost unchanged for both spin channels. Fig.|4]thus 
confirms that edge functionalization always promotes lift- 
ing of the spin degeneracy, and indicates that such effect 
decays slowly as a function of M. 

A remarkable exception are the MS N0 2 -ZGNRs. The 
MS-N0 2 -2x2Z and MS-N0 2 - 4x2Z are half semicon- 
ducting, and the a and (3 gaps are respectively reduced 
to ~90% ~15% of the original gap, while the MS-N0 2 - 
6x2Z ribbon is metallic for both spin channels. 
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NH f 12x3A 



Nedge-12x3A Bedge-12x3A Nbulk-12x3A Bbulk-12x3A 

DOS arb. u. 



6Nbulk-12x3A 
2Pyr-2vac 



6Nbulk-12x6A 
2Pyr-2vac 



Figure 5: (Color online) Spin density of states of a (red/gray solid line) and (3 (blue/black solid line) for functionalized AGNRs 
in comparison with clean ribbon (black dotted line), (a) 12x3 AGNR functionalized in one edge with NH2. (b) N, (c) B edge 
substitutional atoms on the 12x3 AGNR. N and B bulk substitution for the 12x3 AGNR is shown in (d) and (e), where a 
semiconductor-metal transition occurs. (f,g) two Pyridine-like substitutional doping in the bulk of 2x3a and 12x6a presenting 
metal behavior. Arrows indicate impurity levels due to substitutional atoms. The DOS of the clean ribbon is shifted for the 
bulk N and B substitution such that the edge of the valence band of clean and doped ribbons coincide 



E. Atomic substitution 

In semiconducting CNTs, substitutional B (69l.l70| and 
N (70l . [7H act as acceptors or donors, inserting levels 
in the gap or causing a semiconductor-metal transition. 
B substitution on the edge of meta-stable ferromagnetic 
ZGNRs was shown to induce metal-semiconductor tran- 
sitions [62]. It is thus interesting to study in more detail 
the effects of atomic substitutions in GNRs. 

We start by considering O substitution on the edge of 
ZGNRs. By using a 6x2Z ribbon, we find that in the op- 
timized geometry the two C-0 distances are 1.41A. The 
resulting DOS is reported in Fig. [2jg), showing that the 
appearance of new states in the conduction band pro- 
motes the onset of a semiconductor-metal transition. 

To investigate the role of the functionalization density, 
we repeat the same calculation for a 6x6Z ribbon. In 
this case, no semiconductor-metal transition is observed. 
Instead, the spin degeneracy is lifted, and the band gap 
for the two spin channels are 0.449eV and 0.399eV, which 



represent a change of -22% and -31% with respect to the 
pristine ribbon. Overall, the effects of functionaliza- 
tion (Fig. EJf)) and edge substitution (Fig. EJg)) differ 
only for the the position of the extra levels. Since in the 
center of the ribbon the spin density is much lower than 
the edges, the effect of bulk substitutions on the break- 
age of the spin degeneracy is expected to be weaker than 
the corresponding edge substitution or edge functional- 
ization. Fig. [2][h) shows the effect of N substitution in 
the center of a 6 x 5Z ribbon, which corresponds 1.7 at%, 
i.e. IN atom in a lattice of 59 C atoms. This gives a 
semiconductor-metal transition, with a weak spin degen- 
eracy lifting. However, as shown in Fig. EJi) , the intro- 
duction on the edge of the same ribbon of a NH 2 radical 
leads to a much more pronounced difference between the 
two spin channels DOS. In this case, the N substitutional 
atom is placed in the opposite carbon sub-lattice with re- 
spect to NH 2 , and the separation between N and the NH 2 
radical is maximized. This arrangement is chosen with 
the aim to alter the AF state of both sub-lattices. 

We now consider low concentration atomic substitution 
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on the edge and bulk of AGNRs. For edge substitution 
we use the 12x3 AGNR (1 substitutional atom in 71 C 
atoms, 1.4 at%). In this case N impurity states appear 
~1.6 eV below the valence band (Fig. GDJb)) and the B 
states ^0.8 eV above the conduction band (Fig. (He)). 
The impurities modify the gap by +8.7 meV (+1.5%) for 
N and -32 meV (-5%) for B. In the relaxed structure the 
two N-C nearest neighbor distances are 1.33 A and 1.36 
A and the B-C bonds are 1.40 A and 1.43 A. 

Bulk substitution is then investigated in the 12x3 AG- 
NRs. Here we observe that both N and B give impurity 
levels localized around the substituted atoms. This shifts 
the Fermi energy ~0.5 eV inside the former conduction 
and valence band, resulting in a semiconductor-metal 
transition. As indicated by the arrows in Fig. [5jd,e), 
the impurity bands are localized at ~ 0.15 eV from the 
Fermi energy. In the relaxed structure the N-C and B-C 
distances are 1.41 A and 1.49 A, respectively. 

Unlike CNTs [fiH, [7(3, El, substitutional B and N in 
AGNRs do not give electronic states in the band gap, 
even if semiconductor-metal transitions can still be ob- 
served. Further investigations with different substitu- 
tional atoms or geometries are needed to determine if 
and how conventional doping of GNRs can be achieved. 

We finally consider the effect of pyridine-like doping. 
This consists in the substitution of 4 carbon atoms by 3 
nitrogens and a vacancy, as schematically shown in the di- 
agram of Fig.[Hf,g). In CNTs, high density pyridine-like 
doping was found to introduce donor-like states above 
the Fermi e nerg y and to produce semiconductor-metallic 
transitions [72]]. We start by introducing 2 pyridine- 
like impurities in the 12x3 AGNR, corresponding to 9.4 
at% N. The vacancies are arbitrarily positioned near the 
center of the ribbon at 1 row and 2 columns of atoms 
away from each other. This gives 4 localized levels in 
the valence band near the Fermi energy, which induce 
a Fermi energy shift and cause a semiconductor-metal 
transition. We also consider two pyridine-like impurities 
in the 12x6A ribbon, corresponding to 4.4 at% N. These 
are positioned in the nearest possible sites to the center 
of the ribbon at a inter- vacancy distance in the periodic 
direction of 5 columns of carbon atoms. The arrow in 
Fig. EKg) shows the impurity levels in the valence band of 
the ribbon due to the pyridine-like doping. These levels, 
as in the case of higher impurity density, shift the Fermi 
energy and cause a semiconductor-metal transition. 



V. DISCUSSION 

The fabrication of devices with a new material requires 
the understanding all the factors that influence its elec- 
tronic and spin properties. The geometry and the pres- 
ence of various form of chemical modifications are crucial 
for GNRs. The geometry is defined by length, width, 
and edges chirality, whereas chemical modifications can 
be caused by different edge terminations, and by substi- 
tutional atoms in the ribbons' body. 



In this paper, we examined how all these factors af- 
fect the electronic and spin properties of infinite ribbons. 
Our main result is that a strong spin polarization is main- 
tained in all the simulated zigzag ribbons. In particular, 
our calculations show that single edge functionalization 
leads to half-semiconductors with different band gap for 
each spin and can also result in a spin-polarized half- 
semiconductor or in a semiconductor-metal transition. 
This implies the possibility of using single edge func- 
tionalized ribbons to realize spin filtering devices, even 
though single-edge functionalisation may be experimen- 
tally very challenging. On the other hand, the weak de- 
pendence of spin degeneracy lifting and half-metallicity 
on width is encouraging for potential applications, since 
wider ribbons are easier to make. 

We also observe that high concentrations of oxygens on 
the ribbons edges are likely to produce metallic ZGNRs, 
while for lower oxygen concentrations, ZGNRs are ex- 
pected to be half semiconducting. Thus, oxygen is one 
of the most effective elements to promote the closure of 
the electronic gap for both spin channels. Since ribbons 
can be cut from bulk graphene by means of an oxygen 
plasma [1,[78|, this could affect their final electronic prop- 
erties. However, it is likely that OH may also result from 
this treatment. But OH does not close the gap. Given 
the large gaps so far reported [1, [zl], it is thus unlikely 
that these GNRs are just terminated. 

It is known that the electronic structure of AGNRs 
does not present spin-polarization effects [H, [3§] and 
that their gap scales inversely with the width [35|, l39| . 
Our calculations indicate that the gap in AGNRs is ex- 
tremely stable with respect to the edges chemical termi- 
nation. This is due to the absence of edge states, which 
is, on the contrary, the main characteristic of zigzag rib- 
bons. Armchair ribbons seem thus ideal for devices. 
However, our calculations indicate that the presence of 
B and N does not give defect states that pin the Fermi 
energy in the gap, showing that the behavior of AGNRs 
with respect to atomic substitutions is markedly different 
from that of CNTs. Indeed, even if it is well established 
that replacing C atoms with B (N) in semiconducting 
CNTs results in electronic states that pin the Fermi en- 
ergy close to the conduction (valence) band, producing a 
n(p) -doping (fH, [7(3, ED, III E3, [H] , our simulations show 
that in AGNRs the same substitutions lead to completely 
different effects. If the substitution occurs on the ribbon 
edges, the impurity levels are always very far from the 
Fermi energy, and the electronic properties of the ribbon 
are substantially unchanged. If it occurs in the center, 
the defect states appear close to the bottom (top) of the 
conduction (valence) band of the pristine ribbon, but the 
Fermi energy is shifted in the conduction (valence) band, 
triggering a semiconductor-metal transition, rather than 
doping, unlike in CNTs and bulk semiconductors. 

The physical reason for the different effects of chem- 
ical doping in GNRs and CNTs is likely to derive from 
the different boundary conditions to which the electronic 
wave-functions have to obey in the two materials. How- 
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ever, further investigations are needed to fully address 
this point. The reason why edge substitutional atoms 
do not result in a shift of the Fermi energy is related to 
the relative position of defect states and the bandgap of 
the pristine ribbon. Indeed, N edge atoms are attached 
to the ribbon with very localized chemical bonds in an 
sp 3 configuration. This promotes the formation of a deep 
state in the valence band, which is occupied by the ex- 
tra electron carried by the N. On the contrary, N in the 
bulk of the ribbon is in an sp 2 configuration. This gives 
delocalized states in the conduction band. Since no new 
states are formed in the valence band, the extra elec- 
trons coming from N have to be allocated in what was 
the conduction band of the pristine ribbon, with a con- 
sequent shift of the Fermi energy. A perfectly specular 
consideration can be done for B doping. 

Finally, we stress that here we focussed only on 
chemical modification of the ribbons, without taking 
into account additional effects, such as chemisorption 
or physisorption of molecules on the ribbons'surfaces. 
Since adsorbates ca n eas i ly shift the fermi le vel in large 
graphene samples |l08l . 1 1141 . Ill5l . Illd Ill7| . the same 
could also happen for ribbons. This deserves further in- 
vestigation, and may help understand some experimental 
observations, such as the tendency of small ribbons to ap- 
pear positively charged [13, [zll, [ZJ] or the dependence of 
their minimum conductivity on the environment 



ferromagnetic phase of the spin on the edges of zigzag 
nanoribbons is altered by any single edge functionaliza- 
tion, producing half-semiconducting structures. An ap- 
propriate functionalization with the correct radical could 
produce half-metallicity, since different radicals alter the 
anti-ferromagnetic arrangement in different ways. Oxy- 
gen edge substitution and functionalization at high den- 
sity in zigzag nanoribbons produce semiconductor-metal 
transitions. Double edge functionalization of zigzag rib- 
bons can reduce the semiconducting gap, with an anti- 
symmetric arrangement of the radicals energetically fa- 
vored. Nitrogen bulk substitution at 1.7 at% in zigzag 
ribbons also makes them metallic. The gap of armchair 
ribbons is robust against edge functionalization, and N,B 
edge substitutions. However N, B and pyridine-like bulk 
substitution causes semiconductor-metal transitions. 
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